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The concept of magnetic polaron is analysed and developed to elucidate the nature of itinerant 
charge carrier states in magnetic semiconductors and similar complex magnetic materials. By con- 
trasting the scattering and bound states of carriers within the s — d exchange model, the nature 
of bound states at finite temperatures is clarified. The free magnetic polaron at certain conditions 
is realized as a bound state of the carrier (electron or hole) with the spin wave. Quite generally, 
a self-consistent theory of a magnetic polaron is formulated within a nonperturbative many-body 
approach, the Irreducible Green Functions (IGF) method which is used to describe the quasiparticle 
many-body dynamics at finite temperatures. Within the above many-body approach we elaborate a 
self-consistent picture of dynamic behavior of two interacting subsystems, the localized spins and the 
itinerant charge carriers. In particular, we show that the relevant generalized mean fields emerges 
naturally within our formalism. At the same time, the correct separation of elastic scattering cor- 
rections permits one to consider the damping effects (inelastic scattering corrections) in the unified 
and coherent fashion. The damping of magnetic polaron state, which is quite different from the 
damping of the scattering states, finds a natural interpretation within the present self-consistent 
scheme. 
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I. INTRODUCTION 

The properties of itinerant charge carriers in complex magnetic materials are at the present time of much inter- 
est. The magnetic polaron problem is of particular interest because one can study how a magnetic ion subsystem 
influences electronic properties of complex magnetic materials. Recently, semiconducting ferro- and antiferromagnetic 
compounds have been studied very extensiveljiS*^*^. Substances which we refer to as magnetic semiconductors, 
occupy an intermediate position between magnetic metals and magnetic dielectrics. Magnetic semiconductors are 
characterized by the existence of two well defined subsystems, the system of magnetic moments which are localized 
at lattice sites, and a band of itinerant or conduction carriers (conduction electrons or holes). Typical examples are 
the Eu-chalcogenides, where the local moments arise from 4f electrons of the Eu ion, and the spinell chalcogenides 
containing Cr^^ as a magnetic ion. There is experimental evidence of a substantial mutual influence of spin and 
charge subsystems in these compounds. This is possible due to the sp — d{f) exchange interaction of the localized 
spins and itinerant charge carriers. An itinerant carrier perturbs the magnetic lattice and is perturbed by the spin 
waves. It was shown that the effects of the sp — d or s — f exchang 

06/18,9,10^ as well as the sp — d{ f) hybridization , 
the electron-phonon interaction and disorder effects contributed to essential physics of these compounds and various 
anomalous properties are found. In these phenomena, the itinerant charge carriers play an important role and many 
of these anomalous properties may be attributed to the sp — d{f) exchange interaction^"- As a result, an electron 
travelling through a ferromagnetic crystal will in general couple to the magnetic subsystem. From the quantum 
mechanics point of view this means that the wave function of the electron would depend not only upon the electron 
coordinate but upon the state of the spin system as well. Recently, further attempts have been made to study and 
exploit carriers which are exchange-coupled to the localized spinsiiiiii^iiSiii. The effect of carriers on the magnetic 
ordering temperature is now found to be very strong in diluted magnetic semiconductors (DMS) ^i^^ . Diluted mag- 
netic semiconductors are mixed crystals in which magnetic ions (usually Mn^'^) are incorporated in a substitutional 
position of the host ( typically a II- VI or III-V ) crystal lattice. The diluted magnetic semiconductors offer a unique 
possibility for a gradual change of the magnitude and sign of exchange interaction by means of technological control 
of carrier concentration and band parameters. 

It was Kasuya2ii24i24i2iS£ who first clarified that the s-f interaction works differently in magnetic semiconductors and 
in metals^. The effects of the sp-d(f) exchange on the ferromagnetic state of a magnetic semiconductor were dis- 
cussed iQ22i2i2ii2S42&. It was shown that the effects of the sp — d{f) exchange interaction are of a more variety in 
the magnetic semiconductorsiS than in the magnetic metals^iiSLS^, because in the former there are more parameters 
which can change over wide rangesiSiiSiiSiSS. The state of itinerant charge carriers may be greatly modified due to the 
scattering on the localized spinsSSiSSi^iiS^. Interaction with the subsystem of localized spins leads to renormalization 
of bare states and the scattering and bound state regimes may occuri^i^. Along with with the scattering states, an 
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additional dressing effect due to the sp-d(f) exchange interaction can exist in some of these materials^iiSSi^. To some 
extent, the interaction of an itinerant carrier in a ferromagnet with spin waves is analogous to the polaron problem 
in polar crystals^l if we can consider the electron and spin waves to be separate subsystem^SSi. Note, however, that 
the magnetic polaron differs from the ordinary polaron in a few important pointsSSiSSi^SiSi^. 

To describe this situation a careful analysis of the state of itinerant carriers in complex magnetic materials^ is highly 
desirable. For this aim a few model approaches have been proposed. A basic model is a combined spin-fermion model 
( SFM ) which includes interacting spin and charge subsystemsSSiSMii^Si^. 

The problem of adequate physical description of itinerant carriers ( including a self-trapped state) within various types 
of generalized spin-fermion models has intensively been studied during the last decade322ii242ii4i^. The dynamic in- 
teraction of an itinerant electron with the spin-wave system in a magnet has been studied by many authorsSiifliSSiS^, 
including the effects of external fields^. It was shown within the perturbation theory that the state of an itinerant 
charge carrier is renormalized due to the spin disorder scattering. The second order perturbation treatment leads to 
the lifetime of conduction electron and explains qualitatively the anomalous temperature dependence of the electrical 
resistivitj*SS'2i2liiLi&. The polaron formation in the concentrated systems leads to giant magneto-resistive effects in 
the Eu chalcogenidesi*^. 

The concept of magnetic polaron in the magnetic material was discussed and analysed in RefsiSSiSSiSSiiiiiS. The future 
development of this concept was stimulated by many experimental results and observations on magnetic semiconduc- 
^Q].g3j50j5ij52j53_ paramagnetic polaron in magnetic semiconductors was studied by Kasuya^^, who argued on the 
basis of thermodynamics, that once electron is trapped into the spin cluster, the spin alignment within the spin cluster 
increases and thus the potential to trap an electron increases. The bound states around impurity ions of opposite 
charge and self-trapped carriers were discussed by de Gennea^. Emin26, defined the self-trapped state and formulated 
that 

" the unit comprising the self-trapped carrier and the associated atomic deformation pattern is referred to 
as a polaron, with the adjective small or large denoting whether the spatial extent of the wave function of 
the self-trapped carrier is small or large compared with the dimensions of a unit cell" . 

In paperaiSi^^iSi^SiS, a set of self-consistent equations for the self-trapped ( magnetic polaron ) state was derived and 
it was shown that the paramagnetic polaron appeared discontinuously with decreasing temperature. These studies 
were carried out for wide band materials and the thermodynamic arguments were mainly usedS in order to determine 
a stable configuration. Some specific points of spin-polaron and exiton magnetic polaron were discussed further in 
papersamSM^. 

Properties of the magnetic polaron in a paramagnetic semiconductor were studied by Yanase^'^, Kubler^* and by 
Auslender and Katsnelson~i. The later authors§S developed a detailed theory of the states of itinerant charge carriers 
in ferromagnetic semiconductors in the spin-wave (low temperature) region within the framework of variational ap- 
proach. The effect of the electron-phonon interaction on the self-trapped magnetic polaron state was investigated by 
Umehar a^^i^^ . A theory for self-trapped magnetic polaron in ferromagnetic semiconductor with a narrow band was 
formulated by Takeda and Kasuya^. The electron crystallization in antiferromagnetic semiconductors was studied 
by Umehar£t£2i2i, and a dense magnetic polaron state was conjectured to describe the physics of Gd^-xvSi. 
A model of the bound magnetic polaron, i.e., electron trapped on impurity or by vacancjsS was developed in pa- 
pgj.g72j73 rpj^g general thermodynamic model of the bound magnetic polaron and its stability was considered in 
j^gfgj74j75j76^ Then similar models were studied by several author s2i2L22,. 

The state of a conduction electron in a ferromagnetic crystal ( magnetic polaron ) was investigated by Richmond^S 
who deduced an expression for one-electron Green function. Shastry and Matti^ presented a detailed analysis of 
the one-electron Green function at zero temperature. They constructed an exact Green function for a single electron 
in a ferromagnetic semiconductor and highlighted the crucial differences between bound- and scattering-state contri- 
butions to the electron spectral weight. A finite temperature self-consistent theory of magnetic polaron within the 
Green functions approach was developed inpSi. 

Recently, new interest in the problem of magnetic polaron was stimulated by the studies of magnetic and transport 
properties of the low-density carrier ferromagnets, diluted magnetic semiconductors (DMS)^*^i224Si24. The concept of 
the magnetic polaron, the self-trapped state of a carrier and spin wave, attracts increasing attention because of the 
anomalous magnetic, transport, and optical properties of DMS^iSi^Si^LSiSS and the perovskite manganites^2i2S*2i. 
For example, a two-component phenomenological model, describing polaron formation in colossal magnetoresistive 
compounds, has been devised recentlyS^. The paper— includes a detailed analysis of the polaron-polaron interaction 
effects in DMS. 

The purpose of the present work is to elucidate further the nature of itinerant carrier states in magnetic semiconduc- 
tors and similar complex magnetic materials. An added motivation for performing new consideration and a careful 
analysis of the magnetic polaron problem arise from the circumstance that the various new materials were fabricated 
and tested, and a lot of new experimental facts were accumulated. This paper deals with the effects of the local 
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exchange due to interaction of carrier spins with the ionic spins or the sp — dov s — f exchange interaction on the state 
of itinerant charge carriers. We develop in some detail a many-body approach to the calculation of the quasiparticle 
energy spectra of itinerant carriers so as to understand their quasiparticle many-body dynamics. The concept of 
the magnetic polaron is reconsidered and developed and the scattering and bound states are thoroughly analysed. 
In the previous papers, we set up the formalism of the method of Irreducible Green Functions (IGF)^'^. This IGF 
method allows one to describe quasiparticle spectra with damping for many-particle systems on a lattice with complex 
spectra and a strong correlation in a very general and natural way. This scheme differs from the traditional method 
of decoupling of an infinite chain of equationsSI and permits a construction of the relevant dynamic solutions in a 
self-consistent way at the level of the Dyson equation without decoupling the chain of equations of motion for the 
GFs. 

In this paper, we apply the IGF formalism to consider quasiparticle spectra of charge carriers for the lattice spin- 
fermion model consisting of two interacting subsystems. The concepts of magnetic polaron and the scattering and 
bound states are analysed and developed in some detail. We consider thoroughly a self-consistent calculation of quasi- 
particle energy spectra of the itinerant carriers. We are particularly interested in how the scattering state appears 
differently from the bound state in magnetic semiconductor. 

The key problem of most of this work is the formation of magnetic polaron under various conditions on the pa- 
rameters of the spin-fermion model. It is the purpose of this paper to explore more fully the effects of the sp-d(f) 
exchange interaction on the state of itinerant charge carriers in magnetic semiconductors and similar complex magnetic 
materials. 



II. THE SPIN-FERMION MODEL 

The concept of the sp — d{ or d— f) model plays an important role in the quantum theory of magnetismSiifliiSiiSiii. 
In this section, we consider the generalized sp — d model which describes the localized 3(i(4/)-spins interacting with 
s(p)-like conduction (itinerant) electrons ( or holes ) and takes into consideration the electron-electron interaction. 
The total Hamiltonian of the model ( for simplicity we shall call it the s — d model) is given by 

H^Hs + Hs-d + Hd (1) 
The Hamiltonian of band electrons (or holes) is given by 

Hs = J2Yl ^^Aa(^o<y n,„n,^„ (2) 

This is the Hubbard model. We adopt the notation 

k k 

In the case of a pure semiconductor, at low temperatures the conduction electron band is empty and the Coulomb 
term U is therefore not so important. A partial occupation of the band leads to an increase in the role of the Coulomb 
correlation. It is clear that we treat conduction electrons as s-electrons in the Wannier representation. In doped DMS 
the carrier system is the valence band p-holes. 
The band energy of Bloch electrons is defined as follows: 

Uj = N-^ ej:exp[ik{R, - Rj)], 

k 

where N is the number of lattice sites. For the tight-binding electrons in a cubic lattice we use the standard expression 
for the dispersion 

£g — 2 t(aa) cos(fcaa) (3) 

a 

where Sa denotes the lattice vectors in a simple lattice with the inversion centre. 

The term Hs-d describes the interaction of the total 3d(4f)-spin with the spin density of the itinerant carriersii 

Hs-d = -2 ^ = ^ ^[5l>L«fe+9-- + z,S%alak+ga] (4) 

i kq cr 



4 



where sign factor is given by 



and 



z„ = (+or-) jor o- = (t or |) 



n-a _\S-q if = + 

" 1 51,, \in = - 

For example, in DMS^ the local exchange coupling resulted from the -p — d hybridization between the Mn d levels and 
the p valence band / ~ ^^-a ■ magnetic semiconductors this interaction can lead to the formation of the bound 
electron-magnon ( polaron-like ) bound state due to the effective attraction of the electron and magnon in the case 
of antiferromagnetic coupling ( J < ). 
For the subsystem of localized spins we have 



Here we use the notation 



2^ ' ^ 2 

u q 



1 .± 2 



Jy - 5] CXp[zfc(i?, - 



This term describes a direct exchange interaction between the localized 3d (4f) magnetic moments at the lattice 
sites i and j. In the DMS system this interaction is rather small. The ferromagnetic interaction between the local 
moments is mediated by the real itinerant carriers in the valence band of the host semiconductor material. The carrier 
polarization produces the RKKY exchange interaction of local momentsS'S'iSiSiiSi 

Hrkky = — ^ KijSiSj (6) 

We emphasize that Kij ^ \P\ ^ V^_^. To explain this, let us remind that the microscopic modelii, which contains 
basic physics, is the Anderson-Kondo modelii 

H = Y^Y^ tija\^aj„ - V Y^^{at„dj^ + h.c.) 

ij a ij a 

-E,J2Y.^t + lj2Untnt. (7) 

i (T icr 

For the symmetric case U = 2Ed and for U ^ V Eq.© can be mapped onto the Kondo lattice modeU^ ( KLM ) 



Here I ^ The KLM may be viewed as the low-energy sector of the initial model Eq. ( [7| ) . 

We follow the previous treatments and take as our model Hamiltonian expression ^ . For the sake of brevity we omit 
in this paper the U-term (low-concentration limit) . This U-term can be included into consideration straightforwardly 
(see Ref4^). As stated above, the model will represent an assembly of itinerant charge carriers in a periodic atomic 
lattice. The carriers are represented by quantized Fermi operators. The lattice sites are occupied by the localized 
spins. Thus, this model can really be called the spin-fermion model. 
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III. OUTLINE OF THE IGF METHOD 



In this section, we discuss the main ideas of the IGF approach that allows one to describe completely quasiparticle 
spectra with damping in a very natural way. 

We reformulated the two-time GF method^S; to the form which is especially adjusted to correlated fermion systems 
on a lattice and systems with complex spectra. A very important concept of the whole method is the Generalized 
Mean Fields (GMFs), as it was formulated in^'^. These GMFs have a complicated structure for a strongly correlated 
case and complex spectra, and are not reduced to the functional of mean densities of the electrons or spins when one 
calculates excitation spectra at finite temperatures. 
To clarify the foregoing, let us consider a retarded GF of the formSi 

G"^ A{t),A^t') »= -i0{t-t') < [A{t)A^{t')]n = ± (9) 

As an introduction to the concept of IGFs, let us describe the main ideas of this approach in a symbolic and simplified 
form. To calculate the retarded GF G(t — t'), let us write down the equation of motion for it 

ujG{uj) [A, > + << [A, H]^ I A^ »^ (10) 

Here we use the notation << A{t),A^{t') >> for the time-dependent GF and << A \ A''' >>cj for its Fourier 
transform24. The notation [A, B] ^ refers to commutation and anticommutation depending on the value of 77 = ±. 
The essence of the method is as follows^'^: 

It is based on the notion of the "IRREDUCIBLE" parts of GFs (or the irreducible parts of the operators, A and A^ , 
out of which the GF is constructed) in terms of which it is possible, without recourse to a truncation of the hierarchy 
of equations for the GFs, to write down the exact Dyson equation and to obtain an exact analytic representation for 
the self-energy operator. By definition, we introduce the irreducible part (ir) of the GF 

('"^ « [A,H]-\A^ »=« [A,H]^-zA\A^ » (11) 
The unknown constant z is defined by the condition (or constraint) 

<[[A,H]^^\A%,>=0 (12) 

which is an analogue of the orthogonality condition in the Mori formalism^^-^^. Let us emphasize that due to the 
complete equivalence of the definition of the irreducible parts for the GFs ((*'"'<< [A, i/]_|v4^^ >>) and operators 
((''■) [^,i/]_) = ([A, iJ] _)(*'■) we win use both the notation freely ( « A\B » is the same as « (A)(*'')|B » 
) . A choice one notation over another is determined by the brevity and clarity of notation only. From the condition 
()12|l one can find 

_ <P,g]-,At]„> _ Afi 

<[AA%> A/o ^'''> 

Here Mq and Mi are the zeroth and first order moments of the spectral density. Therefore, the irreducible GFs are 
defined so that they cannot be reduced to the lower-order ones by any kind of decoupling. It is worth noting that 
the term "irreducible" in a group theory means a representation of a symmetry operation that cannot be expressed in 
terms of lower dimensional representations. Irreducible (or connected ) correlation functions are known in statistical 
mechanics. In the diagrammatic approach, the irreducible vertices are defined as graphs that do not contain inner 
parts connected by the G'^-line. With the aid of the definition Hll() these concepts are expressed in terms of retarded 
and advanced GFs. The procedure extracts all relevant (for the problem under consideration) mean-field contributions 
and puts them into the generalized mean-field GF which is defined here as 

(w - z) 

To calculate the IGF << [A, H]_{t), A^t') » in ((Tn|l . we have to write the equation of motion for it after 

differentiation with respect to the second time variable t' . The condition of orthogonality (|12|l removes the inhomo- 
geneous term from this equation and is a very crucial point of the whole approach. If one introduces the irreducible 
part for the right-hand side operator, as discussed above for the "left" operator, the equation of motion l|l()(l can be 
exactly rewritten in the following form: 

G = G" + G"PG° (15) 
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The scattering operator P is given by 

P= (Mo)-i( « [A,H]^\[A\H]^ »(''-))(Mo)-i (16) 

The structure of equation ( enables us to determine the self-energy operator M by analogy with the diagram 
technique 

P = M + MG^P (17) 

We use here the notation M for self-energy ( mass operator in quantum field theory ). From the definition H17|l it 
follows that the self-energy operator M is defined as a proper (in the diagrammatic language, "connected") part of the 
scattering operator M = (P)^. As a result, we obtain the exact Dyson equation for the thermodynamic double-time 
Green functions 

G^G° + G°MG (18) 

The difference between P and M can be regarded as two different solutions of two integral equations (|15|l and H18|l . 
However, from the Dyson equation H18|l only the full GF is seen to be expressed as a formal solution of the form 

G= [(G°)-i - Af]-i (19) 

Equation (|19|) can be regarded as an alternative form of the Dyson equation (|18|) and the definition of M provides 
that the generalized mean-field GF is specified. On the contrary , for the scattering operator P, instead of the 
property G^G~^ + G^M = 1, one has the property 

(G°)-i - G-i = PG"G-i 

Thus, the very functional form of the formal solution H19(l precisely determines the difference between P and M . 
Thus, by introducing irreducible parts of GF (or irreducible parts of the operators, out of which the GF is constructed) 
the equation of motion (|l()|l for the GF can exactly be ( but using the orthogonality constraint 112|l l transformed into 
the Dyson equation for the double-time thermal GF (|18|l . This result is very remarkable because the traditional form 
of the GF method does not include this point. Notice that all quantities thus considered are exact. Approximations 
can be generated not by truncating the set of coupled equations of motions but by a specific approximation of the 
functional form of the mass operator M within a self-consistent scheme expressing M in terms of the initial GF 

M w F[G] 

. Different approximations are relevant to different physical situations. 

The projection operator technique has essentially the same philosophy. But with using the constraint (|12|l in our 
approach we emphasize the fundamental and central role of the Dyson equation for calculation of single-particle 
properties of many-body systems. The problem of reducing the whole hierarchy of equations involving higher-order 
GFs by a coupled nonlinear set of integro-differential equations connecting the single-particle GF to the self-energy 
operator is rather nontrivial. A characteristic feature of these equations is that besides the single-particle GF they 
involve also higher-order GF. The irreducible counterparts of the GFs, vertex functions, serve to identify correctly 
the self-energy as 

M ^ G(7^ - G-^ 

The integral form of the Dyson equation p8|l gives M the physical meaning of a nonlocal and energy-dependent 
effective single-particle potential. This meaning can be verified for the exact self-energy using the diagrammatic 
expansion for the causal GF. 

It is important to note that for the retarded and advanced GFs, the notion of the proper part M — {P)p is symbolic in 
nature^^. In a certain sense, it is possible to say that it is defined here by analogy with the irreducible many-particle 
T-matrix. Furthermore, by analogy with the diagrammatic technique, we can also introduce the proper part defined 
as a solution to the integral equation 1)17(1 . These analogues allow us to better understand the formal structure of the 
Dyson equation for the double-time thermal GF, but only in a symbolic form . However, because of the identical form 
of the equations for GFs for all three types ( advanced, retarded, and causal ), we can convert our calculations to 
causal GF at each stage of calculations and, thereby, confirm the substantiated nature of definition (|17l) ! We therefore 
should speak of an analogy of the Dyson equation. Hereafter, we drop this stipulating, since it does not cause any 
misunderstanding. In a sense, the IGF method is a variant of the Gram-Schmidt orthogonalization procedure. 
It should be emphasized that the scheme presented above gives just a general idea of the IGF method. A more exact 
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explanation why one should not introduce the approximation already in P, instead of having to work out M, is given 
below when working out the application of the method to specific problems. 

The general philosophy of the IGF method is in the separation and identification of elastic scattering effects and 
inelastic ones. This latter point is quite often underestimated, and both effects are mixed. However, as far as the 
right definition of quasi-particle damping is concerned, the separation of elastic and inelastic scattering processes is 
believed to be crucially important for many-body systems with complicated spectra and strong interaction. 
From a technical point of view, the elastic GMF renormalizations can exhibit quite a nontrivial structure. To obtain 
this structure correctly, one should construct the full GF from the complete algebra of relevant operators and develop 
a special projection procedure for higher-order GFs, in accordance with a given algebra. Then a natural question 
arises how to select the relevant set of operators {^1,^2, ■■■An} describing the "relevant degrees of freedom". The 
above consideration suggests an intuitive and heuristic way to the suitable procedure as arising from an infinite chain 
of equations of motion (|10|) . Let us consider the column 



where 



A,^A, A2 = [A,H], A3^[[A,HiHi^^.An^[l^^[A,H]:^H] 



Then the most general possible Green function can be expressed as a matrix 

1(4 4 ... Ai)» 



G=« 



/A,\ 
' A2 * 



\AnJ 



This generalized Green function describes the one-, two-, and n-particle dynamics. The equation of motion for it 
includes, as a particular case, the Dyson equation for single-particle Green function, and the Bethe-Salpeter equation 
which is the equation of motion for the two-particle Green function and which is an analogue of the Dyson equation, 
etc . The corresponding reduced equations should be extracted from the equation of motion for the generalized GF 
with the aid of special techniques such as the projection method and similar techniques. This must be a final goal 
towards a real understanding of the true many-body dynamics. At this point, it is worthwhile to underline that the 
above discussion is a heuristic scheme only, but not a straightforward recipe. The specific method of introducing the 
IGFs depends on the form of operators An, the type of the Hamiltonian, and conditions of the problem. 
Here a sketchy form of the IGF method is presented. The aim is to introduce the general scheme and to lay the 
groundwork for generalizations. We demonstrated inM that the IGF method is a powerful tool for describing the 
quasiparticle excitation spectra, allowing a deeper understanding of elastic and inelastic quasiparticle scattering 
effects and the corresponding aspects of damping and finite lifetimes. In the present context, it provides a clear link 
between the equation-of-motion approach and the diagrammatic methods due to derivation of the Dyson equation. 
Moreover, due to the fact that it allows the approximate treatment of the self-energy effects on a final stage, it yields 
a systematic way of the construction of approximate solutions. 



IV. CHARGE AND SPIN DEGREES OF FREEDOM 



Our attention will be focused on the quasiparticle many-body dynamics of the s-d model. To describe self- 
consistently the charge dynamics of the s — d model, one should take into account the full algebra of relevant operators 
of the suitable "charge modes" which are appropriate when the goal is to describe self-consistently the quasi-particle 
spectra of two interacting subsystems. 

The simplest case is to consider a situation when a single electron is injected into an otherwise perfectly pure and 
insulating magnetic semiconductor. The behavior of charge carriers can be divided into two distinct limits based on 
interrelation between the band width W and the exchange interaction I: 

|2/S'| > W] \2IS\ < W 
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Exact solution for the s-d model is known only in the strong-coupling limit, where the band width is small compared to 
the exchange interaction. This case can be considered as a starting point for the description of narrow band materials. 
The case of intermediate coupling, when |2/S'| ~ W, makes serious difficulties. 

To understand how the itinerant charge carriers behave in a wide range of values of model parameters, consider the 
equations of motion for the charge and spin variables. 

[ak,j,Hs]- = efeflfeo. (20) 

[ak.,H,_d]_ = J2iSZqag+k-. + z^Sl^a.+ka) (21) 

q 

= -IN-' ^[25,^_,a;^ap+,x - S+_^ial^a,+,^ - ^^a^+^x)] (22) 

pq 

[SZk,Hs-d]- = -IN''' Y^[2S^_gal^ap+q^ - 5'~_^(a^^ap+qt - aj^^Op+^x)] (23) 

pq 

[SI,H.-d]- = -IN-' E(^^««Pi«P+?T - S^^q4^<^P+<il) (24) 

pq 

[S+,H,]^ = ^ J,iS^^S+_^ - 5^,5+) (25) 
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Q 

From Ea. (|20|l - Ea. (|26|l it follows that the localized spin and and itinerant charge variables are coupled. 
We have the following kinds of charge 

Ofecr, '^ka = aicr^fccr 

and spin operators 

(t'I ^^al^ak+qi; ^^al^ak+qi 
q q 

There are additional combined operators 

bka — ag+fe-CT + ZcrSZqaq+ka) 

q 

In the lattice ( Wannier ) representation the operator bka reads 

ha = (S^'^ai^a + ZaS-a^a-) (27) 

It was clearly shown in Keis&^'SS'^^^ that the calculation of the energy of itinerant carriers involves the dynamics 
of the ion spin system. In the approximation of rigid ion sping-^, i.e.: 

the energy shift of electron was estimated as 



2 

qi 



2 4 ^e(fc)-e(fc-g) 
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The dynamic term was estimated as 



4 ^ e{k) - e{k ~ Q) 2N {2tt)^ J ^ e{k) - t{k ~ q) - I< S'- > +Dq^ 

For the case of rare-earth metals, the electron-magnon interaction in the s—d model within a second-order perturbation 
theory was studied by Liu and Davis^i and by Kim22i2S within the Bogoliubov-Tyablikov^'* GFs method. 
To describe self-consistently the charge carrier dynamics of the s—d model within a sophisticated many-body approach, 
one should take into account the full algebra of relevant operators of the suitable " modes" (degrees of freedom) which 
are appropriate when the goal is to describe self-consistently quasiparticle spectra of two interacting subsystems. An 
important question in this context is the self-consistent picture of the quasiparticle many-body dynamics which takes 
into account the complex structure of the spectra due to the interaction of the " modes" . Since our goal is to calculate 
the quasiparticle spectra of the itinerant charge carriers, including bound carrier-spin states, a suitable algebra of the 
relevant operators should be constructed. In principle, the complete algebra of the relevant "modes" should include 
the spin variables too. The most full relevant set of the operators is 

That means that the corresponding relevant GF for interacting charge and spin degrees of freedom should have the 
form 

f << ai<T|a]^, >> « ai„\S^ » « ai„\Sf » « ai„\a'[^,S^ » << ai<^|a]_^,5'/' >> ^ 

« S!\al^, » « S!\S^ » « S!\sf » « S!\al^,S^ » « S!\al_^,sf » 

« Sr-'la]^, » « Sr'^lS^ » « S^'^isf » « S-"ia]^,S'; » « Sf''\al_^,sf » 

<< Sftti^la]^, >> << Sfai^jSj >> << Sfai^lSf » « Sf ai<^|a]^,5'| >> << S!ai^\a''^_^,Sf » 

V« S'r"a«-a|a]„, » « ST'a^-alSI » « Sr''a^-„\sf » « S~''a^^„\a]^,S] » « 5'r"az-<x|a]_^,5j' »J 

(28) 

However, to make the problem more easy tractable, we will consider below the shortest algebra of the relevant 
operators (a^cr, a^^, 6^0-1 ^lo-)- However, this choice requires a separate treatment of the spin dynamics. 
Here we reproduce very briefly the description of the spin dynamics of the s-d model for the sake of self-contained 
formulation. The spin quasiparticle dynamics of the s-d model was considered in detail in papers^Siiiiffi. We consider 
the double-time thermal GF of localized spins^** which is defined as 

g+~ik;t~t')=«S+it),SZkit')»^-iOit-t') < [S+it),SZkit')]- >= 

/ + 00 
dLoexp{-iLut)g+^{k;uj) (29) 
-CXj 

The next step is to write down the equation of motion for the GF. To describe self-consistently the spin dynamics of 
the s — d model, one should take into account the full algebra of relevant operators of the suitable " spin modes" which 
are appropriate when the goal is to describe self-consistently the quasiparticle spectra of two interacting subsystems. 
We used the following generalized matrix GF of the form^SiMiiffi: 



«S+\SZ,» «S+\aZ,» 
«^t\S-k» «<^k\<^-k» 



= g{k;uj) (30) 



Let us consider the equation of motion for the GF Q{k]uj). By differentiation of the GF << {t)\B[t') >> with 
respect to the first time, t, we find 



UJ 



«St\B»^={'^^'"\^^^'>\+ (31) 



pq 

where 

B = 
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Let us introduce by definition irreducible (ir) operators as 

{Sg) = Sg- < Sq > Sgfi-, (al+qaapa) = ap+go-OptT- < ajcr^po- > '^9,0 (32) 
US',) '''SL, - iSL,) '^'S+) - = {{S^g) - (S^,) ^"-5+) - (0, - (33) 

From the condition (|12|l 

< [((^|) - (SL,) - {cf>g - </.fc-,)5+), >^ 

one can find 



2<S^> 



(34) 



iff =< (5,^) -(5,^) - >; X-+=<5-5+> (35) 



Using the definition of the mass operator Eq. I|17l ) the equation of motion, Eq. (|31l ) , can be exactly transformed to 
the Dyson equation, Ea. (|18n 

g = go + QoMG (36) 

Hence, the determination of the full GF G has been reduced to that of Qq and M. The GF matrix Qq in the generalized 
mean field approximation reads 



SisWxS -SiliVxS, 

where 

R = ni + ^2lN^'^xl (38) 

The diagonal matrix elements Qq^ read 
where 



« S+\SZ, ^ , .r,o\.sn... (39) 



^ ^ 2K^^ + K~+ 

"1 - - ^^(-^0 - J.) - N-^/'J2(-^, - j^_,)_2_±_^ _ - n{) (40) 



Here the notation was used 



^0 

^.^^^ (41) 
xg(fc,c.)=jV-^^ ^^^+V^''^ (42) 



'^p.fe = (w + ep - Ep+fe - A/) (43) 
A/ = 2/5'^ 
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9 9 
e((JCr) = - Zff/S'z 

^(n|+nx); 0<n<2 
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We assume then that the local exchange parameter 1 = 0. In this limiting case we have 

<< S+\S^ >>°= — (44) 

u;-SAJo-Jk)~jis:EM'^-J,-k){2K^'' + K;[+) 

The spectrum of quasiparticle excitations of localized spins without damping follows from the poles of the generalized 
mean-field GF ^ 

u;{k) = S.{Jo - Jk) + 2^ Y.^J, - J,-k){2K'/- + K;+) (45) 

' q 

It is seen that due to the correct definition of generalized mean fields we get the result for the localized spin Heisenberg 
subsystem which includes both the simplest spin-wave result and the result of Tyablikov decoupling as limiting cases. 
In the hydrodynamic limit ^ 0, a; ^ it leads to the dispersion law u!{k) = Dk^. 
The exchange integral Jk can be written in the following way: 

Jfc =^exp(-zfci?,)J(|i?,|) (46) 

i 

The expansion in small k gives^^ 



«5+|5Z,»°=-^ (47) 



2S_ 

uj — uj{k) 

;(fc ^ 0) = (^.(Jo - Jk) + ^ Y.^J, - J,-k){2Kl^ + K-+)) ^ Dk^ 



q 

z q 

riq = ^(fci?,)'J(|-R»|)exp(-igi?i) 

i 

It is easy to analyse the quasiparticle spectra of the {s — d) model in the case of nonzero coupling I. The full generalized 
mean field GFs can be rewritten as 

<< S'+|5'" >>°= f48) 

« (Ju (T_i, >> = , , „ ,, r (49) 

l-^e//(^)xg(fc,^) 

Here the notation was used 

The precise significance of this description of spin quasiparticle dynamics appears in the next sections. 



V. CHARGE DYNAMICS OF THE s~d MODEL. SCATTERING REGIME 



In order to discuss the charge quasiparticle dynamics of the s — d model, we can use the whole development in 
Section 3. The concept of a magnetic polaron requires that we should have also precise knowledge about the scattering 
charge states. By contrasting the bound and scattering state regime, the properties of itinerant charge carriers and 
their quasiparticle many-body dynamics can be substantially clarified. 
We consider again the double-time thermal GF of charge operators24 which is defined as 

gka{t-t') =«aka{t),al^{t') »= -i9{t-t') < [aka{t),al^{t')]+ >= 

/ + 00 
du) exp{—iujt)gka{i^) (50) 
-CXD 
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To describe the quasiparticle charge dynamics or dynamics of carriers of the s — d model self-consistently, we should 
consider the equation of motion for the GF g: 

cj « aka\al^ >>w= 1 + Cfc << aka]al^ »u - 
IN-'/^ « {SZ^a,+k-a + z,Slga,+ka)\ai, »^= IN'^'^ « bkMa »- (51) 

Let us introduce by definition irreducible (ir) spin operators as 

(S^) = S^- < > 6,,o = (52) 

By this definition we suppose that there is a long-range magnetic order in the system under consideration with the 
order parameter < >. The irreducible operator for the transversal spin components coincides with the initial 
operator. 

Equivalently, one can write down by definition irreducible GFs: 



( '"■ << S_^ag+k-a\al^ »^) =<< S_^aq+k-a\al^ » 

Then the equation of motion for the GF gkai^) can be exactly transformed to the following form: 

{lu - e{ka)) « akaWl^ >>^ +IN-^/^ « Cka\al„ »= 1 (54) 

Here the notation was used 

Ck. = b'^,^J2 iS-q<^q+k-. + zASl,) '"^a.+k.) (55) 

Following the IGF strategy we should perform the differentiation of the higher-order GFs on the second time t' and 
introduce the irreducible GFs ( operators) for the "right" side . Using this approach the equation of motion, Ea. (|34|l . 
can be exactly transformed into the Dyson equation Eg. 1)18(1 

9ka{^) = glai^) + 9ka{^)Mka{uj)gka{uj) (56) 

where 

gl,{u) =« a^Ma »"= {u; - e{ka))-' (57) 

The mean- field GF Ea. ((57|l contains all the mean-field renormalizations or elastic scattering corrections. The inelastic 
scattering corrections, according to Ea. (|5?^ . are separated to the mass operator Mka{i^)- Here the mass operator has 
the following exact representation (scattering regime): 

MU^) = A4^;™(a;) 

= ^E(('^''^« SZ^q^k+q-.\S!a+^^_^ »(-)^^') -I- (58) 

qs 



To calculate the mass operator Mka{^)i we express the GF in terms of the correlation functions. In order to calcu- 
late the mass operator self-consistently, we shall use approximation of two interacting modes for Af^^™. Then the 
corresponding expression can be written as 

^L-™(-)-^Ey^;3iz^^^(-^'-^) 

(gk+q.-A^2){^Im « S%\S-'' »^J+5fe+^_,(^2)(^/m« (5^^) *n(^-,) (59) 
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where 

Fi{cJi,uj2) = {1 + N{u;i)^ f{uj2)) 

N{uj{k)) = [exp{(3uj{k)) - 1]-^ 

Equations and form a closed self-consistent system of equations for one-fermion GF of the carriers for the 
s-d model in the scattering regime. It clearly shows that the charge quasiparticle dynamics couples intrinsically with 
the spin quasiparticle dynamics in a self-consistent way. 

To find explicit expressions for the mass operator, Ea. H43|) . we choose for the first iteration step in its r.h.s. the 
following trial expressions: 

gkA^) = Situ - eika)) (60) 
—Im « S^\SZ^ »« z42S,)6{lj - zMq)) (61) 

Here Lo{q) is given by the expression Eg . (|45|l . Then we obtain 



_ 2P < > ^ l-/fc-,.T+^(^(g)) 
^3/2 2^^-s{k-q,^)-u;{q) 

This result was written for the low temperature region when one can drop the contributions from the dynamics of 
longitudinal spin GF. The last is essential at high temperatures and in some special cases. The obtained formulas 
generalize the zero-temperature calculations of Davis and Liu^^ and the approach of papers^^iSSiSS. The numerical cal- 
culations of the typical behaviour of the real and imaginary parts of the self-energy in generalized Born approximation 
were carried out in Refsf^fliSS. 



VI. CHARGE QUASIPARTICLE DYNAMICS OF THE s-d MODEL. BOUND STATE REGIME. 

In this section, we further discuss the spectrum of charge carrier excitations in the s — d model and describe bound 
state regime. As previously, consider the double-time thermal GF of charge operators << aka{t),a-icr{t') »■ The 
next step is to write down the equation of motion for the GF g: 

{lu - e{kG)) « aka\al^ »^ +IN-^^'' « Cka\al, »^= 1 (63) 

We also have 

{uj ~ e{ka)) « aka\Cl »^ +IN-'/^ « C'k^lC'l »^= (64) 

It follows from Eas. H63|l and (|64|l that to take into account both the regimes, scattering and bound state, properly, 
we should treat the operators (ika-,ci\^ and Cka,Cl^ on the equal footing. That means that one should consider the 
new relevant operator, a kind of 'spinor' (^^") ("relevant degrees of freedom") to construct a suitable Green function. 
Thus, according to the IGF strategy, to describe the bound state regime properly, contrary to the scattering regime, 
one should consider the generalized matrix GF of the form 



= G{k;u;) (65) 



Equivalently, we can do the calculations in the Wannier representation with the matrix of the form 




a]„ » «a,^|Cj^>> 

a]a» «a„\c%» 



G{ij-Lo) (66) 
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The form of Ea. H66() is more convenient for considering the effects of disorder. 

Let us consider now the equation of motion for the GF G(fc; to). To write down the equation of motion for the Fourier 
transform of the GF G{k; oj), we need auxihary equations of motion for the following GFs of the form 

{uj ~ e{k + q - a)) « SZgak+q-a\al„ >>^= -IN-^/^ « SZgCk+q-a\al„ »^ - 

p 

P 

ZaN ^/^^^Jp << {S_'^p_^^^^Sp - Sp" S^f^p^g^)aq+k~a\al^ »Lj (67) 
p 

To separate the elastic and inelastic scattering corrections, it is convenient to introduce by definition the following set 
of irreducible operators: 

{S_pS_g ) = S_pS_g — < Sg > S^Q^p (68) 

/rr — cr Qz Q — aQZ \ ir ( Q^'^ Q—ac<z \ 

-(p+9) P P '^-iP+Q)> ~ v'^-(p+g)'^P "^P ^-{p+i})) 

(< So > {Spfi - 6p-q) + {(p-p - 0„(p+5))) 5;^) 

This is the standard way of introducing the "irreducible" parts of operators or GFs^^. However, we are interested 
here in describing the bound electron-magnon states correctly. Thus, the definition of the relevant generalized mean 
field is more tricky for this case. It is important to note that before introducing the irreducible parts, Eq.J^SJ), one has 
to extract from the GF << SZgCk+q-crlal.^ » the terms proportional to the initial GF << SZgak+q-crlal.^ »■ 
That means that we should project the higher-order GF onto the initial oneSi. This projection should be performed 
using the spin commutation relations: [SZg, SZp] = z^N''^/'^ SZ^g+py [SZg, Sl^^ = z_^2N-'^/'^S''_^^^^y 
In other words, this procedure introduces effectively the spin-operator ordering rule into the calculations. Roughly 
speaking, we should construct the relevant mean field not for spin or electron alone, but for the complex object, the 
"spin-electron", or for the operator {SZgak+g-a)- This is the crucial point of the whole treatment, which leads to the 
correct definition of the generalized mean field in which the free magnetic polaron will propagate. We have then 

<< SZqSlpau+g+p-aWl^ »uj^ Z^N-^/"^ « SZlg^pfk+q+p-aWl^ >>t^ (69) 

+ « StpS_gak+q+p-a\al^ >>LJ 

« SlpSZgak+q+p-aWl^ »^ = « (Sip) '"'SZgak+q+p-aWk^ >>c^ (70) 

+ < > Spfi « SZgak+q~a\al„ »u 

<< S_q SZpak+q+pa\o.\„ >>oj — << {S ^g S'^pttk+q+pa) ^^Wka >>lj (71) 

+ < SZgSg > Sp^-g « aka\al„ » 

^^-{p+q)^P ^ ^P '^'^-(p+g))^'=+9-'^l'^fco- »iaJ = « ('5'_(^+q)('S'p) ^ Sp ('^'l (p+,) ) '^^ ^fc+g-o- lOfccr >>w (72) 

+ < 5o > ((5p,0 - Sp-q) << SZgttk+q-alal^ »uj 

Finally, by differentiation of the GF << SZgdk+q-ait)! o^k^riO) >> with respect to the first time, t, and using the 
definition of the irreducible parts, Ea. (|69l ) - Ea. H72l V the equation of motion, Eq.JHH), can be exactly transformed 
to the following form: 

{lj + z„Lo{q) -e{k + q- a)) « SZgak+g-aWl^ >>^ +IN-^/^ < SZg >« akaWk^ »c^= 

<< SZpttk+p-aWl^ >>^ + « Ag\al^ »^ (73) 

p 

where 

Ag = -IN ^^^^^{{S_gS_pak+g+pcr) + z_^{Slp) ^""SZyak+g+p-^a} - 
p 
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p 



O-k+q-tT 



It is easy to see that 



' ^ {uj + ZaLj{q) ~ e{k + q ~ a)) 



IN ^'^1 — \ / ^ \, , « S_pak+p-a\al^ »^ +7 — ■ . . \, TT « AI«L »ii75) 

{uj + z„Ljj{q) - e{k + q - a)) ^ ' {uj + ZaLj{q) - e{k + q - a)) ^' ""^ 

After summation with respect to q we find 

E 7 — ^ f \ ^ n , vT « ^9l4<T (76) 

^ (tj + z<^w(g) - e(/c + g - cr)) 

Then Eg. 1)751) can be exactly rewritten in the following form: 

< S~'' S"^ > 

^ ^ '"^ ^ (1 -/Afc^(w))(cj + zaw(q) -e(fc + g-cr))' '^^ 



+ ^ (1 - /Afc.(c^))(u; + z.a;(5) - e{k + q - a)) "^^^ ^^^^ 

where 

Afe. (^) = ^ E r-: — . ^ \, ^ ^ (78) 

N [u; + ZaCJ(q) - e[k + q - a)) 

To write down the equation of motion for the matrix GF G{k;Lu) Ea. (|65l ) it is necessary to return to the operators 
Cka- We find 

rAr-l/2Y^r < S _g > < jStg) "'(5'|) > ^ j 

^\l-IAM){^ + ^^Mq)~e{k + q-a))^ u;-e{k + qa) > « + 
'''' '^'(l-/Afc<,(w))(w + z<,w(g)-e(fc + g-cr)) uj-e{k + qa) ' 

^\l~IAUujmLu + zMq)-e{k + q-<j))^ Lo-eik + qa) i « ^i^-l^ka »^ + 

<< r >> =Vf < s_qS^ > I < i^-q) ""('S'q) > i I 

" ^^(l-/Afc,H)(^ + ^,^(q)-£(fc + g-(T)) ^ LO-sik + qa) 

y-r « AI^L »^ I «-B,|cL»^ i (80) 

'^^{l-IAka{uj)){uj + z^uj{q)-e{k + q-a)) lu - e{k + qa) ' 

where 
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The irreducible operators Eq. H68I ) , Eq. (|69l ) - Eq. H72I ) have been introduced in such a way that the the operators Aq 
and Bq satisfy the conditions 

< [A,«L]+ >=< >=0 



< [Bq, 4J+ >=< [Bq, C'lh >= (82) 



•'kal+ ' 1^1' ^ka 

The equations of motion, Eas. (j79l l and (|8()n can be rewritten in the foUowing form: 



Vj 

' (1 - IAkaiuj)){uj + Zauj{q) -e{k + q- a)) 
(l + IAkaH) «Bq[al^ » 



(1 -/Afe^(w))(w -e(fc + (7cr)) 



} (83) 



where 



v-r < SZ^S- > (l + /Afc.M)<(5i^)"-(g|)-> 

^^{l-IAkAu;)){LO + zMQ)-£ik + q-<j)) {1 - IAkA^)){Lo - e{k + qa)) ^ ^ ' 

>r-. « Aq[Cl »^ (1 + /Afc^(a;)) « Bq[Cl^ »^ 

(1 - IAkA^)){^ + zMq) -e{k + q- a)) (1 - IAua{Lo)){w - e{k + qa)) ^ 



Xka{^) - (1 „ lAkAo^m^ + zMq) -e{k + q- a)) 

(l + /Afc.H)<(^i^)"-(^|)"-> ^ 
{l-IAu„{Lo)){u;-e{k + qa)) ^ 



Here xtai^) plays the role of the generalized "susceptibility" of the spin-electron bound states instead of the electron 
susceptibility Xg(fc,cj) in the scattering-state regime Ea. H40l ) ( sec also Refs.^^'^'^). Analogously, one can write the 
equation for the GF << Cka\Cl^ ». 

Now we are ready to write down the equation of motion for the matrix GF G{k;uj), Ea. H65l ). after differentiation 
with respect to the first time, t. Using the equations of motion JSUl, (|SS|l . (ISH ) and (ESl)j we find 

hG{k;Lo)=i + Y,Hp)D{p;uj) (86) 



where 



n-( ^ f-f^ ^ f87) 



y«B,[al»^ «B,[Cl»J I- 



with the notation 



Lolq = (1 - IAkaiLo))iLo + zMq) -e{k + q- a)) (89) 

^k,q = |^-:#^^(- - ^(fc + 9-)) (90) 

(1 -f lAkaioj)) 

To calculate the higher order GFs D{p;uj) in Ea.f I86|l . we differentiate its r.h.s. with respect to the second-time 
variable (f). After introducing the irreducible parts as discussed above , but this time for the "right" operators, and 
combining both (the first- and second-time differentiated) equations of motion, we get the " exact" ( no approximation 
has been made till now) "scattering" equation 

G{k; Lu) = G"(fc; uj) + G°(fc; oj)PG^{k; uj) (91) 
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Here the generalized mean-field GF was defined as 

G'ik;co)^n-li (92) 
Note that it is possible to arrive at equation f l91|l using the symmetry properties too. We have 

q 

G = /t(r!t)-i + i3t<|t(^)(f^t)-i 

f7L>t = ^{p)P{pq) = G= (Cl-^i) 
p 

G= (il-i/)t + (f2-i/)^/-i$(p)P(TO)$t(q)(/-i)t(f2-i/)t 

pq 

p^i-'{YHp)P{pqW{qm-' (93) 

pq 

We shall now consider the magnetic polaron state in the generalized mean field approximation and estimate the 
binding energy of the magnetic polaron. 

VII. MAGNETIC POLARON IN GENERALIZED MEAN FIELD 

From the definition, Eq. ( I92|l , the generalized mean- field GF matrix reads 

where 

detCl = uj- e(fccr) - I^N~\l„{uj) 
Let us write down explicitly the diagonal matrix elements Gj^ and Gq^ 

« akMa »"= {detCl)-' = {lu- e{k<7) - N-'xLi^))-' (96) 

The corresponding GF for the scattering regime are given by Eg . ( I57|l . As it follows from Eas.f [^7jl and ( I96|l the 
mean- field GF << aka-\o-l.„ >>" in the bound-state regime has a very nontrivial structure which is quite different 
from the scattering-state regime form. This was achieved by a suitable reconstruction of the generalized mean field 
and by a sophisticated redefinition of the relevant irreducible Green functions! We have also 

« Gfc^lGt^ {detCiy^Lj - e{ka)){Lj - s{ka) - deth))^ (97) 

It follows from Eq. ( I96|l that the quasiparticle spectrum of the electron- magnon bound states in the generalized mean 
field renormalization are determined by the equation 

Ek„ - e(fca) + I^N-\l„{Eka) (98) 

The bound polaron-like electron-magnon energy spectrum consists of two branches for any electron spin projection. 
At the so-called "atomic limit" ( when = ) and in the limit fc ^ 0, w — > we obtain the exact analytical 
representation for the single-particle GF of the form 

« a,.\al »"= ^^^-^ + IS)-' + ^^5^(^ - lis + I))-' (99) 
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Here the notation S and Sz — ^ means the spin-value and magnetization, respectively. This result was derived 
previously in paper. 

However, our approach is the closest to the seminal paper of Shastry and Mattis.^'', where the Green function treatment 
of the magnetic polaron problem was formulated for zero temperature. Our generalized mean-field solution is reduced 
exactly to the Shastry-Mattis result if we put in our expression for the spectrum, Eq. ( 1971) , the temperature T — 

« a..|aL »° |t=o = W- e{ka) - 6^,21^ S-^^^f^^Y^ (100) 

We can see that the magnetic polaron states are formed for antiferromagnetic s — d coupling ( / < 0) only when there 
is a lowering of the band of the uncoupled itinerant charge carriers due to the effective attraction of the carrier and 
magnon. 

The derivation of Ea.f I100|l was carried out for arbitrary interrelations between the s-d model parameters. Let us 
consider now the two limiting cases where analytical calculations are possible. 

(i) a wide-band semiconductor ( l/jS* W) 

Eki ~ efc + I — + (101) 

(-/) ^ (£fc_, - efc + 2/(5 - Sz)) < S+SZ^ > 
N ^ {ek-g-ek + 2ISz) 25 

(ii) a narrow-band semiconductor { \I\S :s> W ) 

^,..cn^| 2(5+l)(5 + 5.) 1 (e,^,-efc)<5+5Z,> 

E,, _ I{S + 1) + ^-^—^^^-^—-^e, + - ^ (102) 

In the above formulae the correlation function of the longitudinal spin components K^^ was omitted for the sake of 
simplicity. Here W is the bandwidth in the limit / = 0. 

Let us now consider in more detail the low-tcmperature spin- wave limit in Eas.f [TIJT|l and ( I102|l . In that limit it is 
reasonable to suppose that Sz ~ 5. In the spin-wave approximation we also have 

<S+SZ, >^2S{1 + N{u;{q))) 

Thus, we obtain (c.f. Refs.S2ig&) 

(i) a wide-band semiconductor ( 1/|5 ^ W) 

E'^l^^'^+^S+'-^Y. 1—— + ^ E , ^''-'~'ir,, Ni.iq)) (103) 

N ^ (cfc - tk-q + 2/5) N ^ [ek~q - ek - 2/5) 
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(ii) a narrow-band semiconductor { \I\S :s> W ) 



- + " + WSTTf^ + ^ i2sl 1) "(25 C"^) 

We shall now estimate the binding energy of the magnetic polaron bound state. The binding energy of the magnetic 
polaron is convenient to define as 

SB = Ski - Eki (105) 

This definition is quite natural and takes into account the fact that in the simple Hartree-Fock approximation the 
spin-down band is given by the expression 

£fci = ek+ IS 



Then the binding energy eb behaves according to the formula 



19 



(i) a wide-band semiconductor ( I/IS* ^ W) 



-1) sr^ i^k-q — efe) 



(ii) a narrow-band semiconductor { \I\S ^ W ) 



-0 



7V^(2S-hl) (25-^1) 

where 

{2PS)^ 1 



^Bl 



f 1-1—1/1 

„ {ek-q-ek-2IS) ~ W ^ ^ 



(108) 

The present consideration gives the generahzation of the thermodynamic study of the magnetic polaron. Clearly, local 
magnetic order lowers the state energy of the dressed itinerant carrier, with respect to some conduction or valence band. 
It is obvious that below Tjv of the antiferromagnet the mobility of spin polarons will be less than of bare carriersiSi, 
since they have to drag their polarization cloud along. Experimental evidence for magnetic polarons in concentrated 
magnetic semiconductors came from optical studies of EuTe, an antiferromagnetiS^. Direct measurements of the 
polaron-binding energy were carried out in Ref^S^. 

VIII. QUASIPARTICLE MANY-BODY DYNAMICS AND DAMPING OF QUASIPARTICLE STATES 

A. Green E^inction Picture of Quasiparticles 

An effective way of viewing quasiparticles, quite general and consistent, is via the Green function scheme of many- 
body theoryS^ which we sketch below for completeness. 

At sufficiently low temperatures, a few quasiparticles are excited and, therefore, this dilute quasiparticle gas is nearly 
a noninteracting gas in the sense that the quasiparticles rarely collide. The success of the quasiparticle concept in an 
interacting many-body system is particularly striking because of a great number of various applications. However, 
the range of validity of the quasiparticle approximation, especially for strongly interacting lattice systems, was not 
discussed properly in many cases. In systems like simple metals, quasiparticles constitute long-lived, weakly interacting 
excitations, since their intrinsic decay rate varies as the square of the dispersion law, thereby justifying their use as 
the building blocks for the low-lying excitation spectrum. 

As we have mentioned earlier, to describe a quasiparticle correctly, the Irreducible Green functions method is a very 
suitable and useful tool. 

It is knownSi that the GF is completely determined by the spectral weight function A[uj). To explain this, let us 
remind that the GFs are linear combinations of the time correlation functions 

1 /■+°° 

FAB{t - t') A{t)B{t') — / duj eyi^[iLo{t - i')]^As(w) (109) 

27r J-oo 

FBAit' - t) =< B{t')A{t) >= — / dLuaip[iLo{t' - t)]ABAiLo) (110) 
Here, the Fourier transforms Aab{^) and AsAi^) are of the form 

Aba{lo) = (111) 
Q-'27:Y,GM-PEn)ii;iB,PMlA^p.n)d{Er, - E„, - lo) 

Aab ^exp{-Puj)ABA{-Lo) (112) 
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Expressions and (|112|l are spectral representations of the corresponding time correlation functions. The quan- 

tities Aab and Aba are spectral densities or spectral weight functions. 
It is convenient to define 

Fba{0) =< B{t)A{t) r dtoAiuj) (113) 

F^b(O) =< Ait)B{t) r dLjeMl3^)M^) (114) 

Then, the spectral representations of the Green functions can be expressed in the form 

G'-iio) ^« A\B »l= (115) 

- / [eMP^')~vW) 

G''{uj):^«A\B (116) 

- / [eMPu^') - VW) 

The most important practical consequence of spectral representations for the retarded and advanced GFs is the 
so-called spectral theorem. The spectral theorem can be written as 

< B{t')A{t) (117) 

1 

TT 

1 

TT 



du! exp[ia;(t — t')] [exp{(3uj) — rj] ^IttiGab + *e) 

< A{t)B{t') >^ (118) 
dcj exp(/3w) exp[ja;(t — t')] [exp(/3a;) — rf\^'^ImGAB{'-^ + *e) 



Expressions (|117|l and (|118|l are of fundamental importance. They directly relate the statistical averages with the 
Fourier transforms of the corresponding GFs. The problem of evaluating the latter is thus reduced to finding their 
Fourier transforms providing the practical usefulness of the Green functions techniqu o??i?^ . 

The spectral weight function reflects the microscopic structure of the system under consideration. Its Fourier transform 
origination is then the density of states that can be reached by adding or removing a particle of a given momentum 
and energy. 

Consider a system of interacting fermions as an example. For a noninteracting system, the spectral weight function 
of the single-particle GF Gk{uj) —<< aka-'jol.^ >> has a simple peaked structure 

Ak{uj) ^ S{uj - Ek) 

. For an interacting system, the spectral function Ak{w) has no such a simple peaked structure, but it obeys the 
following conditions: 

Aki(^) > 0; / Akiuj)duj =< [aka,alj+ >= 1 



kai 

Thus, we can see from these expressions that for a noninteracting system, the sum rule is exhausted by a single peak. 
A sharply peaked spectral function for an interacting system means a long-lived single-particle-like excitation. Thus, 
the spectral weight function was established here as a physically significant attribute of GF. The question of what is 
the best way of extracting it from a microscopic theory is the main aim of the present theory. 

The GF for a noninteracting system is Gk{uj) = (w — e^)"^. For a weakly interacting Fermi system, we have 
G'fe(w) = {u — tk — -Mfe('^))~^ where Mk{u)) is the mass operator. Thus, for a weakly interacting system, the (5-function 
for Ak{Lo) is spread into a peak of finite width due to the mass operator. We have 

Mk{uj ± ie) = ReMk{uj) T ImMk{uj) = Afe(w) T rfc(w) 

The single-particle GF can be written in the form 

Gk{uo) ^{uj- [efe + Afc(L^)] ± Tk{uj)}-^ (119) 
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In the weakly interacting case, we can thus find the energies of quasiparticles by looking for the poles of single-particle 

GF nm 

Lo = ek + Afc(w) ± Tk{uj) 

. The dispersion relation of a quasiparticle 

e(fc) = efc + Afe[e(fc)]±rfe[e(fc)] 

and the lifetime l/F^ then reflects the inter-particle interaction. It is easy to see the connection between the width 
of the spectral weight function and decay rate. We can write 

Afc(w) = (expipLj) + l)-i(-i)[Gfc(cj + ie) - Gk{iJ - ie)] - (120) 

(exp(/3w) + l)-i- — . ^^f^2i2 , ^ 

[tj-(efc-f Afc(w))]2+r^(cj) 

In other words, for this case, the corresponding propagator can be written in the form 

Gk{t) « cxp(-ie(fc)t)exp(-rfct) 

This form shows under which conditions, the time-development of an interacting system can be interpreted as the 
propagation of a quasiparticle with a reasonably well-defined energy and a sufficiently long lifetime. To demonstrate 
this, we consider the following conditions: 

Ak[e{k)] « e(fc); Tk[e{k)] « e(fc) 

Then we can write 

G'fc(w) = jx-r^. (121) 

[c.-<fc)][l-^U,(fe)]+zr,[e(A:)] 

where the renormalized energy of excitations is defined by 

e{k) = ek + Ak[e{k)] 

In this case, we have, instead of ( I121|l . 



Ak{Lu) = (122) 



[exp(/3e(fc)) + 1] [1 ^— kfe)] - 



dtu ''^"^ {uj-e{k))^ + T^{k) 
As a result, we find 

Gfc(t) =« akait);al^ »= (123) 
= -ie{t)exp{~ie{k)t) exp(-r(fc)t)[l - 



A widely known strategy to justify this line of reasoning is the perturbation theory. In a strongly interacted system on 
a lattice with complex spectra, the concept of a quasiparticle needs a suitable adaptation and a careful examination. 
It is therefore useful to have a workable and efficient IGF method which, as we have seen, permits one to determine 
and correctly separate the elastic and inelastic scattering renormalizations by a correct definition of the generalized 
mean field and to calculate real quasiparticle spectra, including the damping and lifetime effects. 

B. Damping of Magnetic Polaron State 

We shall now calculate the damping of the magnetic polaron state due to the inelastic scattering effects. To obtain 
the Dyson equation from Ea.f l91|) . we have to use the relation f ll7|l . Thus , we obtain the exact Dyson equation, Eq.( 

G{k; uj) = G°(k; uj) + G"(fc; uj)MkaG{k; oj) (124) 
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The mass operator has the following exact representation: 



/O \ , ^ 

Here the notation was used 

^^ «A,\Al»P «A,\Bl»P «B,\Al»P ^ «B,\Bl»P ^ 
llfc^(w)=2^{ ^ — ^ \ \ ^ \ -— 1 (libj 

^ '^kp'^kq ^kp^kq ^kpUJf^g ihpilkq 

For the single-particle GF of itinerant carriers we have 

« akaWk, »= {[« aka\al, »") ' - Sfe^H}"' (127) 
Here the self-energy operator 'Skai^) was defined as 

^kA^)-'-^ , , , . (128) 

We shall now use the exact representation, Ea.f I128|) . to derive a suitable self-consistent approximate expression for 
the self-energy . Let us consider the GFs appearing in Eg. (' 1126)) . According to the spectral theorem, Eas.f [TT7|l and 
( I118() . it is convenient to write down the GF << ApjAj »p in the following form: 

« Ap\Al — / -[expipLu') + 1] / diexp(iw't) < AUpit) >p (129) 

^'^ J-oo (^-uj' + ie J 

Then we obtain for the correlation function < A]^Ap{t) >p 

< AlApit) ^ < S^Ci^^_,Ck+p-At)SZ;{t) > + < 4+,_,$%_,$_p_.(t)ap+,_,(t) > (130) 



A further insight is gained if we select a suitable relevant "trial" approximation for the correlation function in the 
r.h.s. of (|130|l . In this paper, we show that our formulations based on the IGF method permit one to obtain an explicit 
approximate expression for the mass operator in a self-consistent way. It is clear that a relevant trial approximation 
for the correlation function in H130|l can be chosen in various ways. For example, a reasonable and workable one can 
be the following "mode-mode coupling approximation" which is especially suitable for a description of two coupled 
subsystems. Then the correlation function < A'^^Ap{t) >p could now be written in the approximate form using the 
following decoupling procedure (approximate trial solutions): 

< s^,cl^g_,Ck+p-At)sz;it) >^ S,^p < S^SZ^it) X Cl^,_^Ck+q-At) > (131) 

< al+k-a'^-q-a'^-p-At)ap+k-a{t) >^ Sg,p < $%_,$_,_,(t) >< aj+fe_,ag+fc-a (t) > 

Here the notation was introduced 

PP' 

{Sz^p,^g^ms^,{t)) '"^ - S^,^t){SZ^^^pj{t)) > (132) 

The approximation, Ea.f I13l|l . in the diagrammatic language corresponds to neglect of the vertex correction, i.e., 
the correlation between the propagation of the polaron and the magnetic excitation, and the electron and magnon, 
respectively. This can be performed since we already have in our exact expression ( I130|) the terms proportional to 
P and J^. Taking into account the spectral theorem, Eas.f lll7|) and f I118f) . we obtain from Eas.f ll2t?|) - f I132|) 

« Ap\Al »P^ ^-^S,,p [ f ^^^^^^ Fi(^i,c.2) (133) 



N 'J J LO — LOl — L1J2 
-^Im « SZq\S^ (V^'™ Ck+q-a\Cl^g^^ »a 



+^<l,PT?y2^'^l' ~ -^l-q')^ f f f _P2(wi,tJ2,W3) 

N ^ J J J LU - LUl - UJ2 - 

q 
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« B,\Bl »'^^ -S,.p / / '-^F,{cjuu;2) (134) 



—Im « {Sl^) - »^,) (— /m « C,+,.|Ct^,, ». 



where 



Fi(a;i,W2)-(l + A^K)-/(c^2)) (135) 

F2(tJl,C^2,C^3) = {1 + NiLUl)){l + N{l02) ~ fiiUs)) ^ NiL02)fiL03) = (136) 

(l + A^(cji)) (l + 7V(cj2)) - (l + N{lji) + N{u;2)') /(cja) 

The functions Fi{llJi,uj2), Ea.f ll35|l . and i^2(wi, W2, 1^3), Ea.f ll37|) . represent clearly the inelastic scattering of bosons 
and fermions. For estimation of the damping effects it is reasonably to accept that 

«Ap\Bl»P:^ «Bp|Aj»f~0 (137) 

We have then 

^ ^^^ «Ap\Al»P «Bp\Bl»P 

^^ «A,\A\»P «B,\Bl»P . 

We can see that there are two distinct contributions to the self-energy. Putting together formulae ( I133f) - ( I138|l , we 
arrive at the following formulae for both the contributions 



1 



qq' 



„i 1 x"^/ /" /" /" dujTdujoduJ:^ „ , n 1 / x 

nL = irjY.^-h' - J,-,'f / / / f2 c^i,c^2,^3VxTi 140 

N J J J UJ-UJ1-UJ2-LU3 i^kq) 



Equations (|124|l , p25|) , (|139|l , and p4U|) constitute a closed self-consistent system of equations for the single-electron 
GF of the s — d model in the bound state regime. This system of equations is much more complicated then the 
corresponding system of equations for the scattering states. We can see that to the extent that the spin and fermion 
degrees of freedom can be factorized as in Eq. ( I131|l , the self-energy operator can be expressed in terms of the initial 
GFs self-consistently. It is clear that this representation does not depend on any assumption about the explicit form 
of the spin and fermion GFs in the r.h.s. of Eas.f lTS^ and f ll40|l . 

Let us first consider the so-called "static" limit. The thorough discussion of this approximation was carried out in 
Refits.. We just show below that a more general form of this approximation follows directly from our formulae. The 
contributions of the GFs, Eas.f ll3J|> and f I134|l . are then 



7^ f r du)iduj2 
W'^^j J u^-u:,-U2 



« AMI / / ^/":""^_ i^i(^i,^2) (141) 



— Im « S_g\Sg >>^i j [ — Im « Ck+q-a\Cl^q^^ », 



>^^.-q ^q-q, ^ .~ q^ , 'J J^_^^^^^ 

q' 
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« B,\Bl ^<5,,, < {SD - > / (142) 



07 — 071 



— /m « Cfc+g^|C^^ 

In the limit of low carrier concentration it is possible to drop the Fermi distribution function in Eqs. ( 11331) - i I140|l . 
In principle, we can use, in the r.h.s. of Eqs. (|135|l and p37|) . any workable first iteration-step form of the GF and 
find a solution by iteration ( see Ref.«24 ). It is most convenient to choose, as the first iteration step, the following 
simple one-pole expressions: 

— !-/m << SZpak+q+p-a\Spa\^^^^^_^ »^=< SZpSp > 5{uj + z„uj.p - e{k + q + p - a)), 

=^Im « (Sip) ^^ak+,+p.\iS;) ^-^al+.+p, ».==< {S%) ^^S;) - > 5{lo - e{k + q + p - a)), 

—Im « SZ^\S^ -z„2 < S^> S{u; + z^Uq), 

TT 

— !-/m << ak+q+p-a\al+q+p-a >>c^= ^{'^ -e{k + q- a)) (143) 

Using Eas. (|129|) - p35|l in (|125|) we obtain the self-consistent approximate expression for the self-energy operator ( 
the self-consistency means that we express approximately the self-energy operator in terms of the initial GF, and, in 
principle, one can obtain the required solution by a suitable iteration procedure ) 

^"^^^ ^ (w^ J2 \uj + z„(ujq - ojp) - e{k + q - p - <t) u + z^tUq - e{k + q + p - a) J 



qp ^ 

T^— <■ ,(l+iV(c^'))/-l 



(144) 



/ < SZlSl > < (5^^) -(5g) - > 

Vo7 — 07' + z^LUq — e{k + q + p — a) 07—07' — e(k + q + per) 

Here we write down for brevity the contribution of the s — d interaction to the inelastic scattering only. For the 
spin- wave approximation and low temperatures we get 

T fo7^~(2^^)V ^ N{u;p){l + N{u;q)) 

Using the self-energy 'Skai^^) it is possible to calculate the energy shift Afe(j(o7) = Re'Skai^) and damping Tka{^) = 
—ImJ^krri^) of the itinerant carrier in the bound state regime. As it follows from Ea. (|145|l . the damping of the 
magnetic polaron state arises from the combined processes of absorption and emission of magnons with different 
energies {ujq — ujp). Then the real and imaginary part of self-energy give the effective mass, lifetime and mobility of 
the itinerant charge carriers 

^ ^ ' - ^""'dlik^^^^^'^^-^^ 

77I 1 1 

g - (^)(-); - = /mSfc^(£(fca))le(fca)=£F (147) 



IX. CONCLUSIONS 



In summary, we have presented an analytical approach to treating the charge quasiparticle dynamics of the spin- 
fermion (s — d) model which provides a basis for description of the physical properties of magnetic and diluted 
magnetic semiconductors. We have investigated the mutual influence of the s — d and direct exchange effects on 
interacting systems of itinerant carriers and localized spins. We set out the theory as follows. The workable and self- 
consistent IGF approach to the decoupling problem for the equation-of-motion method for double-time temperature 
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Green functions has been used. The main achievement of this formulation is the derivation of the Dyson equation 
for double-time retarded Green functions instead of causal ones. That formulation permits one to unify convenient 
analytical properties of retarded and advanced GF and the formal solution of the Dyson equation which, in spite of the 
required approximations for the self-energy, provides the correct functional structure of single-particle GF. The main 
advantage of the mathematical formalism is brought out by showing how elastic scattering corrections (generalized 
mean fields) and inelastic scattering effects (damping and finite lifetimes) could be self-consistently incorporated in 
a general and compact manner. This approach gives a workable scheme for definition of relevant generalized mean 
fields written in terms of appropriate correlators. A comparative study of real many-body dynamics of the spin- 
fermion model is important to characterize the true quasiparticle excitations and the role of magnetic correlations. 
It was shown that the charge and magnetic dynamics of the spin-fermion model can be understood in terms of the 
combined dynamics of itinerant carriers, and of localized spins and magnetic correlations of various nature. The two 
other principal distinctive features of our calculation were, first, the use of correct analytic definition of the relevant 
generalized mean fields and, second, the explicit self-consistent calculation of the charge and spin-wave quasiparticle 
spectra and their damping for the two interacting subsystems. This analysis includes the scattering and bound state 
regimes that determine the essential physics. We demonstrated analytically, by contrasting the scattering and bound 
state regime that the damping of magnetic polaron is affected by both the s-d and direct exchange. Thus, the present 
consideration is the most complete analysis of the scattering and bound state quasiparticle spectra of the spin-fermion 
model. As it is seen, this treatment has advantages in comparison with the standard methods of decoupling of higher 
order GFs within the equation-of-motion approach, namely, the following: 

At the mean-field level, the GF one obtains, is richer than that following from the standard procedures. The generalized 
mean fields represent all elastic scattering renormalizations in a compact form. 

The approximations ( the decoupling ) are introduced at a later stage with respect to other methods, i.e., only into 
the rigorously obtained self-energy. 

The physical picture of elastic and inelastic scattering processes in the interacting many-particle systems is clearly 
seen at every stage of calculations, which is not the case with the standard methods of decoupling. 
Many results of the previous works are reproduced mathematically more simply. 

The main advantage of the whole method is the possibility of a self-consistent description of quasiparticle spectra 
and their damping in a unified and coherent fashion. Thus, this picture of an interacting spin-fermion system on a 
lattice is far richer and gives more possibilities for analysis of phenomena which can actually take place. In this sense, 
the approach we suggest produces a more advanced physical picture of the quasiparticle many-body dynamics. We 
have attempted to keep the mathematical complexity within reasonable bounds by restricting the discussion, whenever 
possible, to the minimal necessary formalization. Our main results reveal the fundamental importance of the adequate 
definition of generalized mean fields at finite temperatures which results in a deeper insight into the nature of the 
bound and scattering quasiparticle states of the correlated lattice fermions and spins. The key to understanding of 
the formation of magnetic polaron in magnetic semiconductors lies in the right description of the generalized mean 
fields for coupled spin and charge subsystems. Consequently, it is crucial that the correct functional structure of 
generalized mean fields is calculated in a closed and compact form. The essential new feature of our treatment is that 
it takes account the fact that the charge carrier operators {al.^,aka) should be treated on the equal footing with the 
complex "spin-fermion" operators {Cl^,Cka)- The solution thus obtained agrees with that obtained in the seminal 
paper of Shastry and Mattis^, where an approach limited to zero temperature was used. 

Finally, we wish to emphasize a broader relevance of the results presented here to other complex magnetic materials. 
The detailed consideration of the state of itinerant charge carriers in DMS along this line will be considered separately. 
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